Universality in the transport response of molecular wires physisorbed onto graphene electrodes 
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We analyse the transport response of large molecular wires bridging graphene electrodes, where the molecules 
are physisorbed onto the graphene sheets. In our study, the sheets are pulled away to vary the gap length and 
the relative atomic positions. The molecular wires are also translated in directions parallel and perpendicular to 
the sheets. We show that the energy position of the molecular resonances is universal for a given molecule, in 
the sense that it is independent of the details of the graphene edges, gaps lengths or of the molecule positions. 
Although the junctions discussed here show a large variability in the zero-bias conductance, we suggest that this 
variability can be reduced drastically by suitably tuning the energy position of the HOMO or LUMO resonances 
to low enough energies. 

PACS numbers: 72.80. Vp,85.65.+h,73.63.Rt,71.15.Mb 



Single- or few-molecule electronics was considered for 
some time a plausible alternative technology to silicon, that 
could be deployed to enable reaching the atomic limit 1 ' 2 . 
However, the picture that has emerged along the past few 
years is that single-molecule junctions show a serious lack 
of mechanical robustness together with a large variability in 
the electrical response. These drawbacks are linked directly 
or indirectly to the few-particle nature of the contacts in these 
devices. For example, the lack of conductance reproducibility 
arises from the large phase space of coupling configurations 
with different conductances. Hence, it is impossible to attach 
a precise value of conductance for a given combination cou- 
pling group/metal. Instead, experimental results are presented 
in terms of statistically averaged conductance histograms 3 . 
The lack of reproducibility is being reduced by the use of a 
new generation of contact groups 4 " 6 , but this improvement is 
marginal. 

Graphene electrodes could be used as an alternative to con- 
ventional metallic electrodes to help solving this lack of re- 
producibility. The use of this material brings about sev- 
eral clear advantages. First, the lack of atomic diffusiv- 
ity in graphene makes each gap geometry stable over time. 
Second, graphene's planar nature enables the visualization 
of the junctions, and possibly the control over the actual 
molecule positioning. Third, graphene allows more eas- 
ily for gateable junctions. A number of theoretical arti- 
cles have indeed proposed bridging graphene electrodes with 
molecules, whose anchor groups would be chemisorbed to 
the electrodes' edgesZ~—. An alternative approach was ad- 
vanced in Refs. (12 13), where fullerene-based dumbbell 
wires where physisorbed onto graphene, or molecular wires 
were physisorbed onto carbon nanotubes. By contacting 
the molecules deep inside the sheets, trouble related to the 
size-mismatch between the physical gap and the length of 
the molecule is avoided from the outset 12 . In addition, the 
conductance variability inherent to the tribological nature of 
graphene's edges is avoided^. Recently it was also proven 
experimentally that stable and gateable molecular junctions 
can be formed by depositing molecules on top of a graphene 
nanogap fabricated by electroburning 14 , where the molecules' 



anchor groups were probably contacting the edges. 

In this manuscript we explore the electronic properties and 
the electrical response of molecular wires physisorbed onto 
graphene leads. Physisorption is a gentler attachment to 
graphene than chemisorption. The distortion of the geometry 
and the electronic structure of both molecule and graphene's 
contact area is expected to be only moderate and mostly re- 
lated to screening effects. Specifically, charge transfer be- 
tween molecule and electrodes is strongly suppressed so the 
contact-induced dipoles are expected to be small. In contrast 
to Ref. OH), we analyze molecular wires with planar an- 
chor groups such as porphyrine or phthalocyanine molecules. 
These groups are linked by polyine chains to form dimer and 
trimer wires, an example of which is shown in Fig. 0). 
The rationale behind our choice is that increasing the con- 
tact surface area could enhance the mechanical stability of the 
junctions and provide an effective conductance averaging. In 
addition, porphyrine wires show a small conductance decay 
with wire length when contacted by gold electrodes 15 .. Por- 
phyrine wires can also show interference and spin-filtering 
effects that make them promising functional units 16 . As op- 
posed to dumbbell molecules, the central region in our wires 
is by geometrical constraints placed always at the same height 
measured from the graphene layers, bringing a desirable de- 
gree of isolation to the functional central unit and, therefore, 
of reproducibility, 

Our main result is a remarkable universality in the energy 
position of the molecular levels, which is completely missing 
for junctions comprising noble-metal electrodes, or for edge- 
contacted graphene junctions. Specifically, we find that the 
energy position of the molecular levels does not shift when 
the molecules are displaced, or when the graphene sheets are 
pulled away to vary the physical gap length. These universal- 
ity property suggests that the conductance variability can be 
reduced by appropriate tuning of the energy position of the 
lowest-lying molecular levels. 

We will discuss our results in terms of the junction shown in 
Fig- 0> where a trimer molecule made up of three porphyrin 
units is linked by butadiynes, e.g.: atomic chains made of four 
carbon atoms. We expect that this molecule will be placed 
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FIG. 1 : (Color online) Extended molecule used to calculate the transport properties. Small light-grey (graphene) and large dark-grey (molecule) 
spheres and small white spheres represent carbon and hydrogen atoms, respectively. Medium-size spheres around the center of each porphyrin 
ring represent nitrogen atoms. The system has a total of 446 atoms. 



roughly at the same height above the sheets (~ 3.2 A), re- 
gardless of its position. Hence, we will discuss the changes in 
the transmission coefficients and in the conductance when the 
molecule is displaced across the plane in directions parallel 
and perpendicular to the edges, or when the graphene sheets 
are pulled away or pushed towards each other. The sheets 
shown in Fig. (fl} are terminated in a zigzag configuration, 
and passivated by hydrogen atoms in a 1-1-1-sequence. How- 
ever we expect that edges will have a reduced impact when 
the anchors attach deep inside the sheets. We used the Den- 
sity Functional Theory code SIESTA 17 , which employs norm- 
conserving pseudopotentials to account for the core electrons 
and linear combinations of atomic orbitals to construct the va- 
lence states. In order to better describe physisorption, we used 
a recent implementation 18 of the van der Waals functional of 
Dion et al 19 . We computed the Hamiltonian and overlap ma- 
trix elements, the potential and the density using a real-space 
grid with fineness defined with an equivalent energy cutoff of 
200 Ry. We relaxed quickly the coordinates to forces smaller 
than 0.05 eV/A using a single-^ basis set (SZ). We took then 
these coordinates as seeds for further coordinate relaxation 
using a single-^ polarized (SZP). We repeated the procedure 
now with a double-^ (DZ) and then with a double-^ polar- 
ized (DZP). We first relaxed an isolated porphyrin ring on top 
of graphene and used those coordinates as a starting point to 
construct the junction by joining the carbon chains and the 
central porphyrin ring. Tests performed with one of the junc- 
tions showed that there are not very large differences in the 
positions calculated with different basis sets. The electronic 
structure is however very sensitive to the choice of the basis: 
we found large differences in the transmission and the den- 
sity of states when the basis set improved from SZ to SZP 
and to DZ. There were not however many changes in mov- 
ing from DZ to DZP, which implies that the electronic struc- 
ture is well converged with DZ. We used then this basis set 
to perform all calculations. The junctions analyzed here were 



made periodic along the direction perpendicular (y) to elec- 
tronic transport (z). We kept only the gamma point in our cal- 
culations because the lateral size of the graphene electrodes 
is large enough to have reduced finite-size effects. Each elec- 
trode has 10 principal layers, where each principal layer has 
two columns of carbon atoms. This was enough to accommo- 
date the porphyrin molecules on top of the sheets and ensure 
that the electronic structure converges to the bulk electronic 
structure at the leftmost and rightmost layers. The graphene 
layers were pulled/pushed along the z— axis opposite to each 
other, so that the physical gap length d increased or decreased. 
We chose a reference distance do between layers equal to 9.1 
A. This corresponded to the separation in which the central 
porphyrin covers entirely the gap. We varied the gap between 
d = -7.6 A and d = 11.2 A in steps of 0.2 A. 

The transport properties were computed using the multi- 
scale transport code GOLLUM 20 . The zero-voltage transmis- 
sion coefficients T(E, V = 0) as the electrodes gap length 
changes from do are plotted in Fig. (0, at intervals of 0.8 A. 
The figure shows that the position of the molecular resonances 
remains almost constant. We have found that this behavior is 
universal. We have observed it for the present junction and 
also for a monomer, dimer or trimer made up of porphyrine 
or phthalocyanine planar groups physisorbed onto graphene. 
According to these results it is therefore possible to attach 
definite values of the molecular levels energy position to a 
given molecular wire. This universality could have been an- 
ticipated in fact: physisorption here carries no charge transfer 
between the molecule and the sheets, nor associated dipole 
moments. This is so because both sides of each contact are 
mostly made up of the same chemical species: carbon. Also, 
the 7r — 7r hybridization between molecular orbitals and the 
electrode states is weaker than for the bonds present in most 
noble-metal/single-atom contacts, and does not have a large 
impact on the nature of the molecule orbitals. We note that 
this kind of universality should also be robust when strong 
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FIG. 2: (Color online) Evolution of the transmission T(E) as the 
separation between the graphene layers d increases from -7.6 A to 
+ 10.8 A relative to do, in steps of 0.8 A. Panels (a), (b) (c) and (d) 
correspond to curves in the intervals [-7.6 A, -3.6 A], [-2.8 A, 1.2 A], 
[2.0 A, 6.0 A] and [6.8 A, 10.8 A], respectively. Decreasing width of 
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FIG. 3: (Color online) Evolution of the zero-bias conductance as the 
separation between the graphene layers along the transport direction 
(z) increases from -7.6 A to 1 1.2 A, relative to do. for various volt- 
ages and temperatures. 



correlations are included. Strong correlations must lead to a 
rearrangement of the electronic spectrum as compared to the 
DFT spectrum. Molecular resonances are shifted in energy 
and in some cases spitted into Lower/Upper Hubbard bands 
and/or a Kondo peak. However, even the strongly correlated 
spectrum should have these features always placed at the same 
energy position regardless of the relative disposition between 
molecular wire and graphene. 

Fig. (0 also shows that T(E) changes by orders of mag- 
nitude at energies E placed between two given molecular res- 



onances, and in particular at the HOMO/LUMO gap. As a 
consequence, the low-voltage, low-temperature conductance 
G = GqT(Ef,V — 0) shows strong variability as a func- 
tion of the relative molecule/graphene position. An example 
is shown in the top panel of Fig. 0, where we plot G as a 
function of d. This panel also shows the room temperature 
conductance, computed using the approximate expression 



G 



e 

w 



J dE T(E, V = 0) (f(E + eV/2) - f(E - eV/2)) 



(1) 

The use of the above formula can be rationalized by noting 
that the voltage drop must be strongly localized at the contacts 
because of the physisorbed bond and should not have a strong 
impact on the transmission coefficient. Notice that, according 
to Fig. (0, finite temperatures lead only to a small reduction 
in the variability of G. In contrast, a finite voltage may have 
a much stronger impact on G, provided the energy integration 
window includes the HOMO/LUMO resonance, whichever is 
closest to the Fermi energy. This resonance, whose energy 
position is universal, will determine the conductance. We plot 
in the lower panel of Fig. G as a function of d for several 
voltages. In the present case, the conductance is controlled 
by the HOMO, which is placed at about 0.75 eV below the 
Fermi energy. As a consequence, for voltages of order 1 .5 eV, 
the conductance variability is reduced drastically. A closer 
knowledge can be gained by plotting G in a decimal, instead 
of a logarithmic scale. The conductance then still varies by a 
factor of 4 if the anchor groups interact with graphene's edges. 
We show G for distances d > do in Fig. (0. Notice that the 
conductance is essentially constant when 2.5 A< d < 8 A, 
which corresponds to gap lengths such that none of the three 
porphyrine molecules overlap with the edges. 

It is relevant to argue now that organic synthetic chemistry 
is now developed to the level where molecules can be fabri- 
cated with a tailored energy spectrum, where the position of 
the HOMO and LUMO levels can be chosen a priori. This 
could be done for molecules isolated in vacuum, solution, or 
possibly deposited onto graphite. Longer linker groups than 
those used in the present study may also be synthesized to 
avoid graphene's edges. 

We mention finally that we have computed the binding en- 
ergy, Eb, for chemical detachment. Typically, we find energy 
barriers of the order of 2 eV, in agreement with a previous 
analysis 21 . However, we find energy barriers between adja- 
cent minima for phthalocyanine monomers sliding on top of 
a graphene sheet Eb ~ 40 meV, which are also consistent 
with a previous calculation 22 . We can make a rough estimate 
of the molecule's drift length r ~ y/4~Dt using an Arrhenius 
equation for the diffusion constant 



D = D Q e- EB/kBT 



(2) 



where we take Do ~ 10~ 4 cm 2 /s. We find that the molecule 
will drift quickly away from the junction except for tempera- 
tures below 5 — 10 K. Notice however that we have not taken 
account of graphene edges, which should help pinning the 
molecule to the physical gap. 
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FIG. 4: (Color online) Evolution of the zero-bias conductance as the 
separation between the graphene layers along the transport direction 
(z) increases from 0.0 A to 1 1 .2 A, relative to do , for a voltage of 1 .5 
V and a temperature of 298 K. 

In summary, we find that physisorbing molecular wires 
with large-area planar anchor groups onto graphene sheets 
shows universality properties whereby the molecular reso- 
nances are placed always at the same energy position regard- 
less of the relative molecule/sheet position and of the phys- 



ical gap length. This pinning is further enhanced if the an- 
chor groups attach deep inside the sheets and do not over- 
lap with graphene's edges. This universality property is a 
consequence of the iso-chemical nature of the physisorbed 
species. Our finding is sustained by DFT simulations of sev- 
eral graphene junctions bridged by monomeric, dimeric and 
trimeric molecular wires composed of porphyrine or phthalo- 
cyanine molecules. We have used a van der Waals functional 
that is suited for the problem at hand. We argue that even 
if the energy position of the molecular resonances is strongly 
rearranged by strong correlations, the new spectra should fol- 
low the same universality behavior. We suggest that the above 
property could be used to solve the conductance variability 
problem inherent to noble-metal-based single-molecule junc- 
tions. To do so, synthetic chemistry could be used to design 
molecules with HOMO or LUMO levels close enough to the 
Fermi energy that the low-voltage conductance is controlled 
by then. We finally find binding energies for molecule detach- 
ment of the order of 2 eV. 
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